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Abstract 

Pointing out the difference between tfie Discrete Nonlinear Schrodinger equation with the classical power 
law nonlinearity-for which solutions exist globally, independently of the sign and the degree of the non- 
linearity, the size of the initial data and the dimension of the lattice-we prove either global existence or 
nonexistence in time, for the Discrete Klein-Gordon equation with the same type of nonlinearity (but of 
"blow-up" sign) , under suitable conditions on the initial data, and some times on the dimension of the 
lattice. The results, consider both the conservative and the linearly damped lattice. Similarities and differ- 
ences with the continuous counterparts, are remarked. We also make a short comment, on the existence of 
excitation thresholds, for forced solutions of damped and parametrically driven, Klein-Gordon lattices. 

1 Introduction 

Some of the most important phenomena in solid state physics (localization of electrons in disordered media, 
crystal dislocation), are described by inherently discrete models. A particular example is the lattice system 

'(i>n - e(0n-l - 2(?!)„ -I- (j)n+l) + ™0n + Vcj)^ + F{(f>n) = 0, (1) 

describing the dynamics of a one-dimensional chain of atoms, harmonically coupled to their neighbors through 
a parameter e and subjected to a nonlinear (on-site) potential. The unknown (j)n(t) stands for the displacement 
of the atom n, while the lattice may be infinite (n G Z) or finite (|n| < N). The parameters m > 0, > are 
related to the mass of the atoms and some possible linear damping effects respectively. The system is well 
known as the Discrete Nonlinear Klein-Gordon equation (DKG). 

It is shown in that the consideration of weak linear coupling force between the particles, in the analysis 
of slow temporal variations of the wave envelope, retains in full the discreteness of the primary model. More 
precisely, in the case of a cubic and quartic potential, the above consideration associates DKG, with another 
famous discrete nonlinear differential equation 

i^n + e{i^n-i - 2V'« + ^n+i) + G{^m) = 0, (2) 

G{un) = P\4'n\'^4'n, whcrc (3 stands for an anharmonic parameter. The lattice system jSJ, is known as the 
Discrete Nonlinear Schrodinger equation (DNLS). Both systems ^ and arise in a great variety of phenomena 
in the context of nonlinear physics, ranging from applications in nonlinear optics to applications in biology. The 
analysis of their dynamics has attracted considerable efforts, see for example P 1^ 151 [T^ [T^ 1^ 1^ 1^ 1^ 
and references cited in these articles. 

As a continuation of our previous works on lattice differential equations |18[I19| . here we consider the following 
discrete Klein-Gordon equation DKG, considered in higher dimensional lattices (n = (ni, n2, . . . , njv) G ^^), 

4>n + iy<Pn - iAd(t))n + TO0„ + i^(<?!>„) = 0, n £ , , t > (3) 

'/'n(O) = (/>„,0, Un{0) = (f)„S- (4) 

* Keywords and Phrases: Discrete Klein-Gordon Equation, lattice dynamics, blow-up, global nonexistence, global existence. 
AMS Subject Classification: 37L60,37L15,35Q55,34C11 
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where 

(Ad0)„GZ« = 0(rii-l,n2,...,rijv) + '/'(ru ,"2-1,. ■ ■ ,njv ) ^ ^ 0(ni ,ri2 1) 

— 2A^ '/'(ni,ri2,...,njv) + 0(tii + I,n2 ,. . .,njv ) 

(5) 

Equation 0, could also be viewed as the spatial discretization of the KG partial differential equation, 

(t)u + Ht- ^(t) + m(l3 + F{(j)) = 0, X eWL^. (6) 

We shall consider as a model for the nonlinearity F, the classical, power-law nonlinear interaction, 
(P) F(s) = /3|sp'^s, seR, /3 = ±1, cr>0. 

The nonlinearity (P) with /? = —1, is usually called as a "non-dissipative nonlinearity" or "blow-up" term, 
since is responsible for the appearance of finite time blow-up and global nonexistence of solutions of various 
partial differential equations, including the KG partial differential equation. The case (] ~ 1, is sometimes called 
as the case of the "dissipative" nonlinearity. 

The motivation for the study of the lattice system ©-©-(P) is due to the following reasons: 

(i) Comparison with continuous counter parts. As it is stated in |24| . "analogy between lattices and nonlinear 
partial differential equations is limited", where the authors refer to the global existence for the DKG lattice 
with interaction F(0„) = /3[(0„ — <j)n-iY + {(t>n — 4>n+i)'^] (FPU lattice) for tr > 2, and the dispersive blow-up 
of the long-wavelength limit, the modified KdV equation, 

+ (jjxxx + PV~'^4'x = 0, 

which occurs for large values of a. In 'ST, the authors consider the case where /? > 0, including nonlinearities of 
the form PF{s) — j3sP^^, for p ~ odd > 3. In 2Sj, the authors consider strongly damped lattices and moreover, 
detailed decay estimates are proved, providing energy decay at exponential rate. 

Another important and characteristic example, is provided by the DNLS lattice and its continuous limit, 
the NLS partial differential equation. It is observed numerically [S]-we refer also to ^] for a simple proof- 
that solutions of conservative DNLS systems with nonlinearity (P) exist globally in the natural phase space, 
unconditionally with respect to the sign of the nonlinearity /3, the degree of the nonlinearity cr, and the size 
of the initial data. Sharp decay estimates of small solutions of DNLS, both for the cases /3 ± 1, are rigorously 
proved in the recent work 133! . 

Similar behavior is shown for the the solutions of the weakly damped DNLS systems 12, : In the unforced 
case, it can be easily seen, that solutions of the weakly damped DNLS systems decay exponentially. For the 
forced system, it can be shown that a global compact attractor exists. 

This is a vast difference in comparison with the NLS-partial differential equation fpi llUll^ for the conser- 
vative case, |34| for the weakly damped case), for which solutions may blow-up in finite time. 

(ii) Comparison with other lattice systems. Our study, takes into account, the following reported observations 
and results: (a) the strong relation between DNLS and DKG systems, through the analysis of slow temporal 
variations (b) the global existence of solutions, for the conservative or the damped DKG system for /3 = 1 
|24l 1251 1331 1371 IHHj (c) the global existence results and decay estimates for small solutions of 33 , in the case 
P = ±1, (d) the unconditional global existence of solutions for the conservative or damped DNLS system 

[31111 ESI. 

We intend to report in this work, further possible differences and similarities between the DKG equation 
and the KG partial differential equation, as well as between DKG and DNLS equations, regarding the issue 
of nonexistence of global solutions. We also check the effect of linear dissipation, to examine if it sufhces to 
prevent global nonexistence in time. 

Close to the approach of our previous works, where we have used techniques from the theory of infinite 
dimensional systems, this time we are selecting some representative energy arguments |^ IT^ I2(il I29j . from 
the extensive bibliography concerning the analysis of evolution equations of second order in time. 

The preliminary section 2, is devoted to a brief description of the functional set-up, and some local existence 
results ^1 Ej. In section 3, we consider the conservative DKG equation. We are based on the abstract 
framework of jl5), to derive some results, on global existence and nonexistence, for negative initial energies 
and initial data of definite sign. The transfer in the discrete setting of the method of implies similar 
decay estimates and description of the interval where the solution becomes unbounded (Theorem I3.1|l . Thus, 
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concerning issue (i), the conservative DKG equation, exhibits similar behavior to that of the continuous counter 
part, since the assumptions on the initial data for either global existence and nonexistence, are similar. However, 
concerning the degree of the nonlinearity cr, there is not any dependence with its growth and the dimension of the 
lattice, as it is observed earlier, from the results on local existence (Section 2.2). This simply derived observation- 
as a consequence of elementary properties of the discrete phase space-seems however to be important, if viewed 
in comparison with the corresponding well known restrictions for the KG partial differential equations. The 
case of small data, is also treated. In this case we transfer some of the arguments of "F", for global existence. Let 
us mention that we don't claim the optimality of the derived bounds in Theorem 13.21 fexceot their interesting 
dependence on the parameters cr, m), of the initial data and the global solution. We refer for this issue, to 
|33j . Especially in [331 Theorem 8], an analysis using discrete analogues of Strichartz estimates is carried out, 
providing decay estimates, which has been tested numerically. However, taking into account their dependence 
on the parameters, the bounds derived by the energy argument, could be interesting in a combination with [331 
Theorem 8], in order to check the rate of decay of the global solution obtained. 

A result of global existence for the linearly damped DKG equation {v > 0), is discussed in Section 4. 
Although the transformation ip — e^^*4>, brings the linearly damped equation to an undamped (nonautonomous) 
form making possibly applicable, the global existence results of Section 3 at least in the regime « m, we 
present also an alternative approach, based on the modified potential well method, developed in |2(i) for the 
semilinear damped wave equation. The potential well argument is applicable through the interpolation inequality 
[33 Theorem 4.1, pg.682]. 



where (•, ■)i2 stands for the £^-scalar product. This inequality, which could be seen as a discrete analogue of a 
Sobolev-Gagliardo-Nirenberg inequality in the discrete setting, has been proved of fundamental importance for 
the derivation of excitation thresholds for standing wave solutions of the DNLS lattice 14, ,35-; ■ Its application 
in the potential well method, reveals that the global existence and energy decay of the linearly damped DKG 
lattice, has a dependence under suitable assumptions on the initial data, on the nonlinearity exponent and 
the dimension of the lattice, through the restriction a > jj. Thus the result of Section 4 which is valid for 
multidimensional lattices, is compatible with those of jHSj, connected with the iV-dependent excitation threshold 
results of See also Remark |4. II giving also a summary of generalizations and possible extensions. 

Therefore, concerning issue (ii), the DKG equation exhibits drastically different behavior in comparison 
with DNLS systems, since there exist data, for which solutions cannot exist globally in time, while concerning 
global existence, restrictions on the size of the initial data as well as to the growth of the nonlinearity and the 
dimension of the lattice can appear. At this point, we mention that the fundamental difference between global 
nonexistence and global nonexistence which occurs by blow-up, has been analyzed in 0). 

We remark that the methods are also applicable, for more general examples of nonlinearities than the model 
case (P). These examples, can be generated by the assumptions of PH^HI- Numerical testing of the results, 
could possibly be of some interest. 

We conclude in paragraph 14. II with a short comment on the existence of stationary solutions for DKG, and 
forced solutions for parametrically damped and driven lattices, motivated by [211 Section VI]. The comment 
(for the proofs we refer to EIEI)) is related with the appearance of energy thresholds, for the existence of 
such type of solutions for the DKG equation, |33l I35| . 

2 Preliminaries 

This introductory section is divided in two parts. In the first part, we provide some brief description of the 
phase space, for the infinite system of ordinary differential equations ©-Q). For details, we refer to |19l I18| . 
We also discuss some basic properties of the discrete operators defined. In the second part, we state the result 
on local existence of solutions of ©-Q- 

2.1 Phase space and properties of discrete operators 

We shall consider for some positive integer TV, the sequence spaces denoted by 
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Between £p spaces the following elementary embedding relation |17l \'M)\ holds, 

e^CF, m\iP<M\i. l<q<p<oo. (8) 

This is a contrast with the continuous analogues of L^'(ri)-spaces when 51 C has finite measure, since in this 
case, (IHl is in the opposite direction. The case p = 2, stands for the usual Hilbert space of square-summable 
sequences, endowed with the scalar product 

((/.,7/;)f2= 0„V«, (f>,^ef. (9) 
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Considering the operator A^; defined by ©, we observe that for any (f> € i' 

\\AM%<mmi, (10) 

that is, Ad : £^ ^ i'^ is a continuous operator. Considering now the discrete operators V+ : i'^ ^ defined by 

(V^V'')neZ« = {V'(ni + l,n2,...,n„) ~ V'(ni,n2,...,njv)} 

~t~ {V'(ni, 712 + 1, ■■■,riAf) 



and V- : ^2 ^ i"^ defined by 



and by setting 



+ {V'(ni,n2,...,njv + l) ~ V'(ni,n2,. ..,"«)} ' (H) 

(V~'0)„eZ« = {^(ni-l,n2,....njv) - ''/'(ni,ri2,. ..,«„)} 

~l" {"0(111,712 — ~ '0(ni,n2,...,niv) } 

+ {0( ni,Ti2,---,7ijv — 1) - i'i )}' (12) 



C^K V')neZ" — ■'/'(ni,n2,...,n^_i,7i^ + l,n^+i,...,7ijv) )> (13) 

(^K 0)7162" = 0(711,712,. ..,7l»_i,7l»-l,71„ + l,...,71jv) " 0( 711,712,. ..,71iv )> (14) 

we observe that the operator A^, satisfies the relations, 

N 

(-Ad0i,02)£2 = ^(V+0i,V+^2)£2, for aU 1^1,^2 (f5) 

K=l 

(V+0i,02),2 = (V'i,V-02)£2, for all 7/^2 G^^ (16) 

Relations (|15l) . are sufhcient in order to establish that A^ : i"^ ^ defines a self adjoint A^ < operator 
on Examining the graph norm of the operator A^, we see that D{Ad) = X = i"^ . We shall also use the 
following bilinear form, and induced norm 

JV 



(01,02)^2 := ^(V+0i,V+02)£2 +m(V'i, 02)^2, (17) 

re=l 

N 

■■= Ell^«^lll+HI0ll.'- (18) 



We denote by if the Hilbert space with scalar product (|17|l and norm 1)181) . The usual norm of i'^ and l|18|l are 
equivalent, since 



^2 < 11011,2 < (2iV + m)||V'||,^2. (19) 
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2.2 Local existence of solutions 

As for the DNLS lattice jj^li we shall formulate Q-Q, as an abstract evolution equation in ^|^. In the 

case of the DKG lattice, we may set A :— Ad — m, and check that the operator 



1 



z 



B:=(^^ _^ j, DiB):=^i^l^ j : z,u;eii,Azee'}, (20) 
is a skew-adjoint operator, since 

and £>(B) = £'(B*) = £2 x £2 := £2- The operator B, is the generator of a an isometry group T{t) : R C{(-2)- 
It was shown in |il9i Lemma 2.1], that from the nonlinear function (P), we may define a nonlinear map (still 
denoted by F) F : ^ £'^ for any < a < 00. Setting 





^ = ( i ) , ^{^) 

equation Q can be rewritten as 



^ = BV' + J^(V')- (21) 

For fixed T > and ({0„,o}nez, {0n,i}nez)'^ = (0o>i)'^ := ^0 e £l, a function ^ := (0, </>) e C([0,T],^i) is a 
solution of (|2HI . if and only if 

V'(i) = T(0^o + / T(t - s).F(^(s))ds 5(OV'o, (22) 



Our local existence result can be stated as follows: 

Proposition 2.1 We consider the DKG I3)-0)-(P), with v >{). We assume that ipo G £\. Then, there exists 
a function T* : £\ ^ (0, 00] with the following properties: 

(a) For all e £\, there exists ip G C{[0,T*{-)po)),£^) such that for allO<T < T*{4>o), -0 is the unique solution 
of mj) m C([0,T],^2) i^^^ii posedness). 

(h) For allte [0,r*(V'o)); 

r*(Vo)-^> ^^^^^^^^^^ ■■=Tn, R = 2mt)\\,., (23) 



where L denotes the Lipschitz constant associated with the nonlinear operator J-' : £2 ^ £2- Moreover the 
following alternative holds: (i) T*(iPq) = 00, or (ii) T*{'ipQ) < 00 and 

lim ||'0(t)|L2 =00, (maximality). 

TTT*(Vo) ^ 

(c) T* : £2 {0,00] is lower semicontinuous. In addition, if {ipno}neN is a sequence in £\ such that ■0no V'o 
and if T < r*(V'o); then S{t)ipQn S{t)ipQ in C([0,T],£^) (continuous dependence on initial data). 

It is important to remark further continuity properties of the solution, which will be essential for the 
application of the energy arguments, in the next sections. We have the following 

Corollary 2.1 The solution of I3)-0)-(P), is such that (p G C2([0, T]; £2) n C([0, T]; ^2o-+2) allO<a < 00. 

Proof: It follows from Proposition 12.11 equation Q , the continuity of the operator : £'^ ^ £'^ and the 
Lipschitz continuity of the nonlinear map F : £'^ ^ £^ , that solutions of ©-©-(P), are such that G 
C2([0,r];£2), since 

0n + !^0„ = (Ad0)„-m0„-F(0„) = O, 71 eZ^, ,ie[0,r]. 

From the embedding ©, we have that £^ C £2o-+2 j^j. ^^^y < cr < cx) and [351 Vol. II/A, Proposition 23.2 (c)& 
(h), pg. 407], implies that L''^ ([0, T]; f^) c ^''^([O, r];f2<T+2) f^^. 1 < ra < ri < 00. Therefore 

<?!)eW:={0 : (/)eLP([0,r];^2)^ ^g:^?Q0,T];^2^+2)| \<p,q< 00. 

Then from ES". Vol. II/A, 23.13a, pg. 450], we conclude that e C([0, T]; ^^^+2). o 
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3 The conservative DKG equation 

In this section, we discuss the global in time solvability of the conservative DKG equation (^^ = 0). Our 

analysis is based on the study of the total energy functional 



HW-^ll0Wlll+UEl|V^</.(t)|||,+m||0||||-^-^ E |0„(t) 

U=l J riSZ" 



2a+2 



The functional H, is the Hamiltonian of the infinite lattice of nonlinear oscillators: it is conserved as long as 
the unique solution of (jSJ-lEl exists, that is 

W(0o,0i) :=H(0) = H(i), <e[0,T], foraUO<T<r*. (24) 

The case of non-positive initial energy. We discuss first, the issue of global nonexistence of solutions. We 
will show that the abstract framework developed in ^21) can be adapted for the study of the DKG ©-O-(P)- 
The differentiability of the functional 



U=l J ^ neZN 



2a+2 



(25) 



is needed. This can be done by using and working exactly as in |18l Lemma 2.2]. The proof makes use of 
a discrete analogue of dominated convergence theorem 0. 

Lemma 3.1 Let (p £ f' and < cr < oo. Then the functional l^25\) is a C^{f,R) functional ana 

<v'(0),v^>= - 



' N ^ 
E E (^t'l>)n{^t^)n+m E Mn\ 



E \<t>n\^"Mn^ i^ef. (26) 



We seek for some A > 2, satisfying the key condition of (151 Inequahty (2.6), pg. 455]. Using and ll^ . 
the definition of the equivalent norm if (ITHll . ((T^ and ||SJ), we get that 

(V'(<^),0) - AV(0) - - l) II0III + (l - im'lM, (27) 

Setting A = 2ct + 2, we observe that 

(V'(^),0)-AV(^) >a||0||| >0. (28) 
Using Corollarv l2.1l for the time derivatives of fi{t) = ||0(i)|||2, we have the relations 

^l'{t) = 2((/),,/)),2, 

ti"{t) = 2||0||2, +2(V'(<^),0) >2||0||2, +4(a + l)V(0). (29) 
which can be derived by using Q, the Hamiltonian 7i and 125(1 . We have 

^l'{o) = 2{(t>o,(^>l)p, m(o) = ||0oII,'2. 

We assume nonpositive initial energy 

n{0) < 0. (30) 
Now, from the conservation of energy H24|) and (|29|l . (|30|) . we get 

^^'\t)>2m%+A{a + l)(^^ml~niO)^ >{2a + A)m%- (31) 
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On the other hand we have that 

i^i'it)r<i^,{tfm%. (32) 

The differential inequaUty 



which can be written also in the form 



,"it) >^ 1 ^ (34) 



[m'(OF- 2 /i(i)' 

is derived by combination of 132|l and (|31|l . As in |15[ pg. 455], we may conclude with the following observations: 
Assuming that /^'(O) > 0, by integrating H33() . we obtain the Ricatti inequality 

^Ji'{t) ^i'{Q) 

^ , .„^.^ + 2 > 0, (35) 



with blow-up time estimated by 

2 

On the other hand, assuming that /i'(0) < 0, by integration of (|34|l . we get that 

implying that /i'(t) < 0, for all t > 0. We summarize in the following theorem, whose proof is identical to those 
of [151 Lemma 2.1 & Lemma 2.2]. 

Theorem 3.1 We consider the DKG equation i0)-0'-(P), with (3 — ~\ and < ct < -f-cxo. (i) We assume 
that the initial data are such that 7i(0) < and ((/)o, 0i)f2 < 0. Then the unique solution (j) G C'^{[0,+oo),£'^), 
and satisfies the decay estimate 

\\mu^<\\MU<i+st)-K 

(a) We assume that the initial data are such that Ti.(0) < and (0o, (t>i)i2 > 0. Then the solution blows up (in 
the sense that fi(t) = \\u{t)\\'j2 becomes unbounded) on the finite interval (0,r*), with 

= ^H^. (36) 



Global existence results for sufficiently small initial data. In the present paragraph, we show that 
solutions still exist globally in time, if we assume conditions only on the size of the norm of the initial data. 
This time, we choose to provide a "discrete version" of [51 Proposition 6.3.3]. Inspection of the proof, will 
reveal that discreteness may have some additional interesting effects, in comparison with the continuous case 
(see R.emark |4.1ll . The proof could be viewed as a simple alternative to [231 Theorem 8]. The strength of [331 
Theorem 8], is also on the detailed decay estimates of the solution and its connection with excitation thresholds, 
for the existence of standing waves (see ). However, the conditions needed for the result presented here, seem to 
be consistent with the excitation thresholds of paragraph 14. II since there are also independent of the dimension 
of the lattice and the growth of the nonlinearity exponent. 

Theorem 3.2 We consider the DKG equation 0)-0)-(P), with (3 = —1 and < (t < +oo. We consider also 
the function 

0ix) = Mox'+'' - X, x€[0,+^), Mo = ^-j^— -, (37) 
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and we set —p := Ta\n9{x) < 0. Assume that the initial data ipQ = (0o, 0i) £ £2 are sujficiently small in the 
sense 

m\i ■■= m\% + m\% < mm |i, ^r^Si^^^} ' ^''^ 

Then (f> £ C^([0, +00)), and satisfies the estimate 

sup||0WII|. + ||0WII|<^^:5:^^ (39) 

Proof: We shall consider the function 

Kt)^\\m\i + \mmr (40) 

It is convenient to rewrite the conservation of energy (|24|l . as 

flit) ^ A^(O) - 2J(</>o) + 23m), Hm) = E \Mt)\'''+'- (41) 

By using the embedding relation and the equivalence of norms (|19(l . we easily derive the inequality 



< 



1 



2a+2 



< , .L.^i M(0"+^ (42) 



{a + l)m'^+i " 
1 

Thus, from iglj and (g^l we get that 

fi{t) < m - 2J(0o) + -—^—^ii{t)'^+' . (43) 

The same arguments f (l42|l for 0o) can be applied, in order to show that inequality H43(l . holds for the initial 
value /i(0). We have that 

m - 23iM < m + j^^^^^i^i^r^'- (44) 

A first assumption on the size of the initial data is 

M(0) < 1. (45) 

Then, from it follows that 

m(0)-2J((/)o) <Mim(0), Mi = 1 + Mo. (46) 

The rest, will follow by repeating the arguments of Proposition 6.3.3, pg. 85-86]: The consideration of the 
function 6{x), is motivated by inequality 1431) . It appears that for all a € (0,p), there exist XonUa, such that 
< Xa < Ua, satisfying 9{xa) + a = 9{ya) + a = 0. Moreover it holds that a < Xa < a{l + <j)/<J- Combining 
(|43|l and (|46() . we observe that 

e{fi{t)) + Afi/i(o) > 0, 

for all t S [0,T]. Thus, assuming that 

MMO) < p, (47) 

we may set a = A/i/i(0), in order to apply the aforementioned argument, for the function 9{x) and the point 
Xa- We obtain that 

fi{t) G [0,a;Mi^(o)) U (2;mi/.(o),+oo). (48) 
Now the continuity of the function /z(i), implies that 

M(t) < ^Ah^^io) < i^MiA*(0), (49) 

for all t £ [0, T). Inequalities H45I) and H47() . enforce us to choose /i(0) < min(l, pAI^^), hence (|38ll on the size of 
the norm of the initial data. The last part of the rhs of inequality H49|l . is the estimate l)39|l on the corresponding 
global solution, o 
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4 The linearly damped DKG 

In this section, we consider the Unearly damped DKG equation ©-©-(P)- Let us note that the transformation 
ip e~^*(j>, brings the equation into the nonautonomous undamped DKG equation 

4>n - (Arf0)„ + m^(/)„ + /3F(0„, t) = 0, ne Z^, 

where = m — ^ and F(t, 0„) = e~'"^*|(/)„p°'0„. Theorem 13 . 21 seems to be appHcable (after some modifica- 
tions) at least in the regime m >> i/^, which is a physicaUy interesting case. However, we decide to present a 
global existence and energy decay result, for the equation 

4>n+'y^n - (Arf0)„+/3|0„|2-0„ =0, 71 eZ^, t > 0, (50) 
(/>«(0) = 0,1,0, <i>n{0) ^ (Pn.l- (51) 

based on the the "modified potential well" , introduced in (26j for the study of the linearly damped semilinear 
wave equation. The result makes use of a discrete analogue of the Sobolev-Gagliardo-Nirenberg inequality, 
which was proved in |35j . in the context of existence of excitation thresholds for the DNLS equation with 
power nonlinearity. The potential well arguments have their origin in j28| . The result will demonstrate that 
under certain assumptions on the size of the initial data, the system features global existence and energy decay, 
assuming sufficient growth of the nonlinearity, depending on the dimension. A detailed energy decay estimate, 
is also derived, see 

We recall first, the discrete interpolation inequality, which was proved in 35 . 

Theorem 4.1 Assume that a > j^. Then there exists C* > 0, such that for all cj) £ holds the interpolation 
inequality 

E i^"i'^^'^^*f E I'^-i') Eiiv^<^ii?- (52) 

We rewrite for convenience in the case m — 0, the energy quantities 

1 1 ^ 1 

m = -||0(i)||2,+_^||V+<^(t)||2, -—^ ^ |0„(i)|2-+2^ (53) 

K.= l „GZ« 
1 ^ 1 

Vim = -^T.\\^tm\% + ^^T. \^-(*)\'^'''- (54) 

Multiplying equation H5()|l. by (p in the scalar product, and by using H15|l . we derive the relation of dissipation 
of energy 

n{t)=n{0)-v f \\4>{s)\\j2ds, (55) 

which hold as long as the local solution exists. Moreover, we observe that 

-V{t)<H{t). (56) 
With the help of the inequality (|52() . we shall proceed to the proof of 
Theorem 4.2 Let fi — —1, in in and we assume that the nonlinearity exponent satisfies (J > . Setting 

Po ■■= - {i^WMl + -no) + 2(</'o, 
we assume that the initial data \51\) . satisfy the conditions 

1 ^ 

PO < and E W^Uo I \% >\\M I'-t' ■ (57) 
Then the solution of <j) G C^{[0,+oo),e^), and limf^+oo \\<l){t)\k' = 0- 
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Proof: The modified potential well in our case, will be defined as 

W£2 := the interior of the set |(/) e f\K{(j3) = J2\\Vt4'\\p " > o| ■ (58) 

To the fact that Wp is well defined for ct > , is a consequence of the discrete interpolation inequality (|52|l 
(see also [221 Pg- 329, (1.13)], for the continuum case). 

We argue by contradiction. Let us assume that there exists T* > 0, such that 

(/)(t) e Wf2 for aU t e [0,T*), and (/)(T*) G dWe2. 

Then it follows that K{(I){T*)) = 0, and 0(r*) ^ 0. For every t e [0,T*], it follows from (jsHJ, that 

N ^ N 

K — l At — 1 

Multiplication of (|50|) . by (j> in the scalar product, implies the identity 

2-mj{t)),.-2\\m\i + ,.-^^^^^^^ (60) 

K.— 1 

Since K{(l){t) > 0, for aU < i < T*, we obtain from that 







,, , / M ,9 1 , 



< u\\M\l+2{-\\m\% + -\\m\e-}+2iq^o,Mi2+2 I MsWi^ds, (61) 







Moreover, from (|59|l . we have that 

K— 1 

for all < ^ < T*. Now, from the inequality (|62|1 and the identity (|55|1 . we get the estimate 

^110(^)11'^+^ f ms)\\%ds<n{Q). (63) 







Then, inserting (|63|l into 1)61(1 . we get that 



||0(t)||,'2 < ^|j/||0o||,'2 + ^H(O) + 2(</.o,0i),2| :=p2 (64) 
for aU < t < T*. With the bound at hand, we use the inequahty (j5^ . to obtain 

E ^^^/^oEll^fc-^WII'^' o<i<r*. (65) 

We set 

A* '■— C^cPq ■ 

Hence, by using (|65|l . we derive that for t = T*, 

N 

AT AT 



K—l K—l 

N 

= (l-A,)El|V^0(T*)||,^. >0, 



4 = 1 
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under the assumption A* < 1, justifying the first of the conditions on the initial data described in H57|) . The 
contradiction imphcs that T* — oo. Now, since JC{4){t) > 0, for all t G [0, +oo), the dissipation of energy identity 
(j55|l . implies that 



11^ 1 

:,\\m\% + :,J2\\^tm\\i^ < n{o) + -—mt)\\'i\i 



K— 1 



K=l 

Thus, we have that 

N 



1 ^ 



K=l 



which implies the bound 



^ll<^WIII^ + ^El|V+0W||^. <^^i^^, for all <e[0,+cx)). (66) 

With the aim to estimate the quantity /J K{(j){s))ds, we integrate once again H6()(l with respect to time, to get 
the identity 

Jo Jo 

+ V f {^{s),cj,{s)),.ds^ [ ms)\\ltM. (67) 
Jo Jo 

Therefore, from and (|H7jl . we get 

i^(0(s))ds < / ||0(s)||2,ds + |(0(t),</)(i))£2| + |(0o,</.l),2| 

Jo 

/ ms),c^{s))p\ds. (68) 



By using H19|l . and the bounds H63(l and (|64|l . for the second and the last term of the rhs of (|68|l . we obtain the 
estimates 

< \[\\m\i + \\m\\i}< ^^^^^^\ (69) 

/ 1(0(5)7 0('S))£2|'^s - ^ 1 \\ii-'^)\\iA\(t>{s)\\pds 
Jo Jo 

1/2 

< {i^poH{Q))^/V/\ (70) 



On the other hand, H55|l implies that 

d 



^^{Hmm+t)}<nm))- (71) 



We shall integrate H71|l with respect to time, and we will insert the relation 
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(obtained from H58|l and (|25|l 'l. to get the inequality 
{l + t)H{t) < H{Q)+ f n{s)ds 

2 Jo 







< no) + ^ I ms)\\%ds- I v{Hs))d. 



s 







1 / ,, : , M,9 , 1 / w , cr 



< nio) 

^ Jo 

- ^^^^^Uo + X A*(2a + 2) X ^W^))^^- (72) 

where A* = 1 - A* > 0. Finally, from {T^J, the estimates lIHSl) and ^^-^7^, we obtain the energy decay 

estimate 

Hit) < L,^+L,^- (73) 

The constants Li, L2 are found to be 

A*+0' I^ + 3\,, A* +(7/ On,, ,,\ 

L. = ^^^(^PoH(0))V^. (75) 

Let us note e proof remains valid for the case m > 0, under similar assumptions on the initial data, under 
slight modifications, o 

Remark 4.1 (Nonlinearity exponents). A review of the results in Section 3, implies that unique solutions 
exist locally (and globally under suitable conditions on the initial data), independently of the growth of the 
nonlinearity and the dimension of the lattice. The same holds for the conditions on global nonexistence of 
solutions. This observation, although a simple consequence of (jHl, the equivalence of norms of l"^ and of 
the "auxiliary" space of "discrete derivatives" i\ (a discrete analogue of the Sobolev space H^{MJ^)), looks 
interesting if one compares the corresponding results for the KG-partial differential equation © with (P). 
We recall a representative result, stated in |26l Section I, pg. 327], 9 . It is well known that local existence 
depends on the growth of the nonlinearity and the dimension of the domain, since unique local solutions exist 
if < cr < l/(iV — 2) and < cr < 00 only for the cases = 1,2, a consequence of the Sobolev embedding 
theorems. Similar restrictions hold for the case of a bounded domain of M^. 

However, in the global existence result of Section 4, for the linearly damped lattice, it appears a dependence 
not only on the size of the initial data but also on the nonlinearity exponent and the dimension of the lattice. 
This is not only in contrast with the DNLS lattice, but it also demonstrates an interesting difference with the 
continuum DKG partial differential equation: while in the continuum case the nonlinearity exponent is required 
to be smaller than a critical value in the case A > 3 (cr < 1/(A^ ~ 2)), in the discrete case and for initial data 
satisfying l|57|) . the nonlinearity exponent has to he greater than a critical wa^ite-depending on the dimension of 
the lattice-cr > 2/A^. Thus, for global existence and energy decay of some initial data in the linearly damped 
DKG lattice, sufficient growth of the nonlinearity is required. Considering for example the equilibrium solutions 
of the linearly damped DKG, and taking into account the results of jSHl Theorems 8 & 9, pg. 1855] for the 
ID-lattice (see also "37, Remark, pg. 1852]), the result of Theorem 14.21 seems to be also a verification of a 
Weinstein type result on higher dimensional lattices for the DKG, on the instability of nontrivial equilibrium 
solutions (and existence of excitation thresholds) for this type of solutions. 

(Generalizations). The arguments described in Sections 3 and 4, could be repeated for more general nonlin- 
earities than the model case (P). We refer to those of Proposition 6.3.3, Proposition 6.4.1] and those with 
the structure assumed in |15l I29| . For further differences with the continuous problem, we note as an exam- 
ple, that the proof of Theorem 13.21 indicates that assumptions on the growth and the sign of the function 
G{x) = F{s)ds, as those of 9, Proposition 6.3.3 & Remark 6.3.4] for small, do not appear in the discrete 
setting, for global existence of small solutions. 

The results hold also for DKG equations, assuming Dirichlet boundary conditions. However, the equivalence 
of norms in the finite dimensional spaces, allows for a transfer of the calculations of 0], to the finite dimensional 
DKG equation, and to produce a finite time blow-up result (rather than global nonexistence of solutions). 
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The assumptions on generalized discrete operators, which are not necessarily discretization of the Laplacian, 
described in pg- 347-348], and which give rise to the their self-adjointness properties, and definition of 
appropriate equivalent norms |3H1 pg- 349], imply that the results of the previous sections, can be extended to 
DKG equations, involving such operators. 

Other, possibly interesting improvements, could consider (a) the case of nonlinear damping. For this case, 
the references ^| 1201 123 123 j should be valuable tools, (b) Discretizations of reaction-diffusion equations [711^. 
involving "blow-up" type nonlinearities. 

4.1 Forced solutions for parametrically damped and driven lattices and equilibria 
of the DKG equation. 

The short discussion is motivated by |21l Section VI], which considers the driven damped chain, of non. linear 
oscillators 

m(l)n - e((/',i-i - 20„ + (j)n+i) + mul sin0„ + vcfin = / cos{2uJet) sin0„, n G Z. (76) 

Expanding sin(/>„ in the Taylor series, one obtains a cubic DKG equation. The constant loq, is the frequency 
of small amplitude vibrations in a well of a substrate potential, while uje denotes the frequency of the driving 
force /. When the force / and the damping v are small, one may try to look for forced solutions of the form 

(/)„(t) = cos(a;e< + B), n^Z. (77) 

Under the approximation described in |21l pg. 3169], the resulting equation for the real function ■!/;„ describing 
the wave envelope, reads as 

-e(?/>„_i-2V'„H-V„+i) + f^Vn = AV'^, 71 gZ, (78) 
where the constants fJ, A, are given by 

n = ^to(w^-cj^) 

A = sin(2e) 

and 8 denotes the phase shift. 

The simple proofs of the results of [191 118| . obviously can be applied for the rigorous justification of the 
existence and nonexistence of nontrivial forced solutions (equation H78|l . is exactly the same with the one satisfied 
by the complex function ipn, in the ansatz for the breather solution e'^*?A„ of the DNLS equation). The same 
is valid, for the problem of nontrivial equilibria, of the DKG equation ©-(©-(P), satisfying the equation 



-2e--/cos(e) 

27We 



Theorem 4.3 A.(Dirichlet Boundary Conditions- Finite lattice) Let e > 0. For any > there exists a 
nontrivial solution of ^76]) . 

Al. Consider the DKC equation (0)-|^-(P). For any m > 0, there exists a nontrivial equilibrium solution. 
B. (Dirichlet Boundary conditions or infinite lattice) Let e > 0. There exist no nontrivial forced solution for 
\7b)j , of power less than 

P„in(r!,A) - ^. (79) 
B.l There exist no nontrivial equilibria for 0)-|^-(P), of power less than 

l/a 



Pmin(™,'T)= ( . (80) 
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The proof, is made by a discrete analogue of the Mountain Pass Theorem It is vahd in the case of a 
finite dimensional problem when A = const ^Hl, and can be extended to infinite dimensional lattices when the 
anharmonic parameter is an element of an appropriate ^^'-space |18| . The extension of the method in order to 
cover the case A — const in infinite lattices possibly needs the concentration compactness arguments of |35| . 

On the other hand, the non-existence results-obtained by a fixed point argument- are valid for the infinite 
dimensional lattices, both for the case of site dependent and site independent anharmonic parameter. Moreover 
it provides a simple proof on the existence of excitation thresholds both in the case of DNLS and DKG cases, 
|35[I33| . Alternatively to those proved in |2Fj-depending on the dimension of the lattice and the exponent of the 
nonlinearity)-the threshold (|79|l . depends on fl or m (or the frequency of the standing wave of DNLS equation) 
and the exponent of the nonlinearity, see jl8l Theorem 2.5]. 

Theorem l4.3B . is valid for any nonlinearity which defines a locally Lipschitz map on i'^, and similar thresholds 
can be obtained. Another interesting example is the so called saturable nonlinearity, F{s) = j^^^2- We refer 
to jTJ for the saturable DNLS problem. 

In conclusion, the assumptions on the initial data needed for the results of Sections 3 & 4, imply the 
instability of the non-trivial equilibria (or standing waves) for the DKG equation ©-©-(P), for such data. 
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